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Kirchhoff’s kinetic analogy relates the deformation of an incompressible elastic rod to the classical 
dynamics of rigid body rotation. We extend the analogy to compressible filaments and find that 
the extension is similar to the introduction of relativistic effects into the dynamical system. The 
extended analogy reveals a surprising symmetry in the deformations of compressible elastica. In 
addition, we use known results for the buckling of compressible elastica to derive the explicit solution 
for the motion of a relativistic nonlinear pendulum. We discuss cases where the extended Kirchhoff 
analogy may be useful for the study of other soft matter systems. 


Analogies to dynamical problems have been used to 
simplify the physics of various condensed-matter systems, 
ranging from the deformation of elastic bodies to the 
order-parameter profile across an interface between co¬ 
existing phases. A particularly well known example is 
Kirchhoff’s kinetic analogy [1]. In this theory the three- 
dimensional (3D) deformation of a slender elastic rod is 
reduced to the bending deformation of an incompress¬ 
ible curve, representing the mid-axis of the rod. This 
problem, in turn, is analogous to the dynamics of a rigid 
body rotating about a fixed point, where the distance 
along the curve and its local curvature are analogous, 
respectively, to time and angular velocity. When the fil¬ 
ament is confined to a two-dimensional (2D) plane (the 
celebrated Euler elastica [2] ), the equation of equilibrium 
coincides with the equation of motion of a physical pen¬ 
dulum [1, 19]. 

In the examples above the elastic system was reduced 
to an indefinitely thin, incompressible body, whose equi¬ 
librium shape follows the trajectory of a classical dynam¬ 
ical system. In the present work we show that relaxing 
the incompressibility constraint introduces terms akin to 
relativistic corrections to classical dynamics. Within this 
analogy, the compression modulus, Y, plays the role of 
the relativistic particle’s rest mass, and the bendability 
parameter (Y/B) 1 ^ 2 = h ~ x , where B is the bending mod¬ 
ulus, is analogous to the speed of light. The limit of an 
incompressible rod (h —> 0) corresponds to the nonrela- 
tivistic limit. 

Despite the relevance to real systems, including com¬ 
pressible fluid membranes [4], there have not been many 
studies of compressible elastica (see [5] and references 
therein). Following these works, we consider the 2D de¬ 
formation of a compressible filament, represented by a 
planar curve of relaxed length L. The same model ap¬ 
plies to thin elastic sheets, as well as fluid membranes 
[6], provided that they are deformed along a single di¬ 
rection. The deformation away from the flat, stress-free 


* ozzoshri@tau.ac.il 
t hdiamant@tau.ac.il 


state is parametrized by the angle <j>(s) and compressive 
strain y(s), as functions of the relaxed arclength s along 
the curve, s G [0, L\; see Fig. 1(a). We denote the com¬ 
pressed arclength by s, such that 7 = ds/ds and the total 

deformed length is L = ds = q( s)ds. 
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FIG. 1. (a) Deformed and undeformed configurations. The 
straight line (solid, purple) represents a relaxed rod of initial 
length L, and the curved line (dashed, gold) represents a de¬ 
formed rod of total length L. Their arclength parameters are 
s and s respectively. The local in-plane and out-of-plane de¬ 
formations are respectively accounted for by the compression 
field 7 = ds/ds and the angle 4>(s). (b) An illustration of the 
discrete model. The initial, zero-energy configuration (solid, 
purple) consists of N rigid bonds of rest length a and zero 
joint angles. A higher-energy state (dashed, gold) is obtained 
by compression of each bond, a —> 7 , 1 a, and/or a change of 
each joint angle, </>„. The continuous limit of this model yields 
the picture presented in panel (a). 

To obtain the energy functional which keeps the bend¬ 
ing and compression contributions independent, it is in¬ 
structive to start from a discrete model (see Fig. 1(b)). 
Consider a chain of N jointly connected compressible 
rigid bonds of rest length a. The chain’s configuration 
is parametrized by a set of bond strains 7n (1 < n < N) 
and a set of joint angles <j> n (1 < n < N — 1). The 
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energy of a given configuration has a compression con¬ 
tribution, E c = (Y/2a) X^ = 1 (7n a a ) 2 i and a bend¬ 
ing contribution, E b = (B/2a) ii^n+1 - </>n) 2 ■ For 
the sake of the analogy we add a potential energy of the 
form, E p = aJ2n=i'y-nV((j)n), where V{cj>) is an angle- 
dependent potential. This choice of potential energy is 
not artificial. For example, if we include an external 
work on the chain, coupling the force P exerted on the 
boundaries with the chain’s projected length, we have 
E p = Pa cos (f) n . i.e. V(<f>) = P cos </>. We now 

take the continuum limit, N oo and a —> 0 such that 
Na —> L and na —» s, getting, 


E [ 0 (s), 7 (s)] = E b + E c + E p 
2 


= f ds 


f (S) +V i - i)2+7V ^ 


(1) 


TABLE I. Mapping between the parameters of compressible 
elastica and relativistic dynamics 


Elasticity 

Relativity 

relaxed arclengtli a , s 

proper time, r 

compressed arclength, s 

laboratory time, t 

tangent angle, <j> 

angle coordinate, <j> 

curvature, re = dip/ds 

angular velocity, w = d<j>/dt 

compression field a , 7 = ds/ds 

Lorentz factor, 7 = dt/dr 

compression modulus, Y 

minus the rest energy, —me 2 

bending modulus, B 

moment of inertia, ml 2 

bendability, h~ x = ^JYjB 

speed of light, c 


a More accurately, the analogue of r is 70s, where 70 is a 
constant prefactor (see text). 


To obtain the equations of equilibrium for the filament 
one should minimize E with respect to y(s) and <f>(s). 
This gives the two equations, 


Y( 1 -1) + V = 0, 


d 2 cj> dV 


= 0 . 


First integration yields, 


( 2 ) 

( 3 ) 


M = ?(S) -yh-i) 2 - 7 V' = co„ 8 t, (4) 


which depends on the boundary conditions at s = 0, L. 

The mathematical analogy to relativistic dynamics is 
revealed once we transform from the relaxed arclength, 
s, to the compressed one, s, and redefine the potential as 
V = —V. Equation (3) then turns into 


b 4 (44j + — 

ds \ ds J d(j) 


( 5 ) 


This equation is analogous to the Euler-Lagrange equa¬ 
tion for the one-dimensional motion of a relativistic par¬ 
ticle [7], provided that y(s) coincides with the Lorentz 
factor. To obtain 7 we use eqn (4) to eliminate V in 
eqn ( 2 ), and transform to s, which gives, 

7 = yuy (6) 

\/l + h 2 {dcj> / ds) 2 

This expression indeed resembles the Lorentz factor up 
to the constant prefactor, 70 = y/l + 2 H/Y. Absorbing 
this prefactor in the potential V makes eqn (5) identical 
to the relativistic one, with the correct Lorentz factor. 
This completes the analogy. The mapping between the 
variables and parameters of these two systems is summa¬ 
rized in Table I. 

We now use the mapping to identify a new symmetry 
of compressible elastica. First, recall that the equation 
of equilibrium of incompressible elastica is invariant to 


the addition of a constant curvature re, making a flat 
configuration cylindrical. The kinetic analogue of this 
symmetry is the Galilean invariance of the equation of 
motion of a classical free particle to a boost by a constant 
velocity. We are after the corresponding symmetry for 
a compressible elastic filament. It is natural to try the 
analogue of a Lorentz boost. To make the connection to 
a free particle we first remove the potential energy from 
eqn (5), V (</>) = 0, and then divide by B"/ 0 . The resulting 
equation of elastic equilibrium, 


d<p 


ds \ + h 2 {d(j)/ds ) 2 ds 


= 0, 


( 7 ) 


is identical to the equation of motion of a free rela¬ 
tivistic particle. With the Lorentz boost in mind, we 
immediately identify the transformation of coordinates 
(s/h,<t>) —> (S/h, 4>), which leaves eqn (7) unchanged: 


(= ( 7 ~ r y flK \ [s/ h \ 

V </> J 7 J \ $ J ’ 


( 8 ) 


where re = const, and 7 = l/\/l + /i 2 re 2 . The transfor¬ 
mation is a rotation of the material coordinates, {s/h, </>), 
by an angle 9 = sin -1 ( 7 /ire), turning a flat, relaxed con¬ 
figuration into a cylindrical, compressed one; see Fig. 2. 
The rotation leaves the line element on the material- 
coordinates plane unchanged, ds 2 /h 2 + d(j) 2 = dS 2 /h 2 + 
d' L 2 . This reparametrization invariance defines a con¬ 
tinuous family of solutions to the equation of equilib¬ 
rium, containing all the cylindrical, compressed configu¬ 
rations of the filament. These configurations satisfy local 
mechanical equilibrium but differ in their total energy. 
(Obviously, the flat relaxed configuration always has the 
lowest energy.) The actual equilibrium configuration for 
a given problem is selected by boundary conditions. To 
our knowledge, this symmetry of compressible filaments 
has not been recognized before. 

An important physical distinction should be made, 
however, between the two sides of the analogy. In the 
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FIG. 2. Rigid rotation of the coordinate system (s/h,rf)) —► 
( S/h , <f>) by an angle 9. A flat elastic filament of length L is 
described by the blue line. As viewed from the new frame, 
the total length is compressed by a factor 7 , and the angle 
coordinate, <E>, is linearly changing, resulting in a cylindrical 
compressed shape of radius 1/k and length 7 L. 


elastic problem there is a unique zero-energy state, and 
one naturally specifies the length L and the boundary 
conditions with respect to that relaxed state; the com¬ 
pressed length L is determined by minimization. By con¬ 
trast, in relativity there is no preferred reference frame, 
and the duration of the experiment and the boundary 
conditions are given in the laboratory frame, not the 
proper one. This difference leads to the extra prefactor 
70 above. Thanks to this prefactor the elastic strain can 
be either compressive or dilative, whereas the relativistic 
Lorentz factor must be dilative . 1 Note that, because of 
the minus sign in the mapping Y —me 2 , time dilation 
corresponds to elastic compression. 

Despite the physical difference, by absorbing the factor 
70 in the potential, the two problems become mathemat¬ 
ically equivalent. Thus, a known solution in one system 
can be readily transferred to the other. We shall now use 
known results for the buckling of compressible elastica to 
write down an explicit solution to the equation of motion 
of a relativistic nonlinear pendulum. (A solution to this 
dynamical system was previously given only in implicit 
parametric form [ 11 ].) 

To do so, we should specialize to the potential V(ft) = 
Pcos</>. Substitution of V in eqn (2) gives 7 = 1 — 
( P/Y) cos (f>. These expressions for V and 7 turn eqn (3) 


1 If we formulated the elastic problem and its boundary conditions, 
unphysically, based on the ultimate deformed length L, the in¬ 
formation on the unique reference state would be lost, and the 
two problems, the elastic and relativistic, would become identi¬ 
cal. This can be readily shown by transforming eqn ( 1 ) to s and 
minimizing with respect to 7. 


d 2 6 P 2 

B-— + Psin</>—— sin</>cos</> = 0, (9) 

ds 2 1 

which is the equation of equilibrium for a compressible 
elastica under a uniaxial force P [5, 8 ]. The exact solution 
to eqn (9) was derived in ref. [5]. For hinged boundary 
conditions, d(j)/ds = 0 at s = 0, L, it is given by 


<f>(s) = 2 sin 1 


(1 - {qkh) 2 / 7 o ) 1/2 cd (y^fcs, m) 

q -—-— 

[l - {qkh .) 2 /70 cd 2 (^ 70 ks, m)} 


( 10 ) 

Hereafter we use the conventional symbols for the var¬ 
ious elliptic functions, cd, a m, K and n, as defined in 
ref. [9]. In addition, k = y/P/B is the wavenumber of 
the buckled linear elastica. We have defined three addi¬ 
tional parameters which depend on the boundary angle, 
</> 0 = </>(s = 0): 


q = sin( 0 o/ 2 ) 

70 = 1 — ( kh) 2 cosc^o 
m = q 2 [l + (qkh) 2 ] /V (11) 

The parameter 70 = 7 (s = 0) gives the boundary strain, 
which coincides with 70 as defined above. The actual 
wavelength (periodicity) A of the nonlinear buckled elas¬ 
tica is 


A = 4K(m)/(^yok). (12) 


The total rod length matches a half wavelength, L = A/2. 
This expression relates </o to the force P and the system 
parameters via eqn (11). The total deformed length is 
given by L = f Q [1 — (kh) 2 cos (j>]ds, yielding, 


L/L = 1 - (kh) 2 


2(1-q 2 ) 


H (q 2 ,m) 
K (to) 



(13) 


It is readily shown that the three known solutions for in¬ 
compressible nonlinear elastica (h —>■ 0 ) [ 1 ], compressible 
linear elastica (cf >0 <C 1) [5], and incompressible linear 
elastica (h — 0 ,^o 1 ) [ 1 ], are obtained from eqn ( 10 ) 

in the respective limits. 

To apply the analogy, we first transform eqn (9) to the 
deformed arclength s using equations (2), (4) and ( 6 ). 
This leads to 


d ( 1 d(j>\ k 2 . 

wv , -rr H -Sin <j) = 0 

ds y y /1 + h 2 (dcf>/ds ) 2 ds J 70 

Next, substituting 

{h 2 ,fc 2 / 7 o }^{-£ 2 /c 2 , w 2 }, 


(14) 


where u> is the natural frequency of the corresponding 
linear pendulum and £. the pendulum length, the equation 
of motion for a relativistic pendulum [ 10 ] readily follows: 

4- ( . 1 + 0 J 2 siiuj) = 0. (15) 

dt \ ^l-(P/c 2 )(dfi/dt) 2 dt j 
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From eqn (10) we then immediately write down the ex¬ 
plicit solution for the pendulum motion in terms of its 
proper time, 


</>(t) = 2 sin 1 


[1 + (quoP/c) 2 } 1 / 2 cd (cot, to) 

^ [1 + (quoi/c) 2 cd 2 (cot, to)] 1 / 2 . 


(16) 


where, from eqn ( 11 ), q = sin(i/> 0 / 2 ) is related to the 
initial angle of the pendulum, and m = q 2 {\ — (u>£/c) 2 (l — 
q 2 )]. Similarly, the proper period of the pendulum is 
obtained by analogy to 7 qA of eqn ( 12 ), 


T T = AK(m)/oo, (17) 


and the laboratory period from L/("/ 0 L) of eqn (13), 


T t /T T = 1 + 2M/c) 2 (1 


I 2 ) 


H (q 2 ,m) 
K(m ) 



(18) 


Finally, we can transform the solution (16) to laboratory 
time through 

t = r + ——— ( colic ) 2 [ll(g 2 , am(wr, to), to) — lot] , 

OJ 

(19) 

In Fig. 3 we show the resulting motion of the pendulum 
in both proper and laboratory times. The dilation of the 
period of the relativistic pendulum (solid curve), com¬ 
pared to the nonrelativistic one (dash-dotted), illustrates 
the compression of the buckling wavelength in the com¬ 
pressible elastica compared to the incompressible case. 



t,T 


We note that a similar factor to the elastic 70 is en¬ 
countered in certain relativistic scenarios [11, 12]. We 
note also that the usage of the original strain field to de¬ 
rive the elastic 7 looks similar in spirit to the “einbein” 
used to simplify the analysis of actions in different phys¬ 
ical contexts [13]. Unlike those analyses, however, 7 here 
is an actual physical strain, the choice of which is not 
free but dictated by elasticity. 

The extended Kirchhoff analogy presented in this 
Communication should be useful for studying various 
soft matter systems. In principle, compressibility cor¬ 
rections are negligible for sufficiently slender (and there¬ 
fore bendable) objects. There are several indications, 
however, that the compressible-elastica model may be 
qualitatively predictive beyond its strict limits of valid¬ 
ity. For example, it produces negative compressibility 
(decrease of the force P with increasing confinement) for 
bendability h of order 1 [5]. Although the model is 
invalid for such values of /i -1 , this result is in qualita¬ 
tive agreement with recent experimental and numerical 
results for buckled beams [14]. In ref. [4] Diggins et al. 
studied the buckling of fluid membranes under uniax¬ 
ial confinement using molecular dynamics simulations. 
To improve the analysis of the simulation data beyond 
the incompressible Euler elastica, they used, on purely 
empirical grounds, a model which is very similar to the 
one presented here, obtaining good fits to the numerical 
results. The reason for the unexpected success of the 
compressible-elastica model in these two examples is un¬ 
clear at present. Compressibility is essential in the case 
of two-dimensional deformations with non-zero Gaussian 
curvature. The present ID model, as far as we know, has 
no applicability in such cases. It is relevant, however, 
to deformations with cylindrical [15, 16] or conical [17] 
symmetries. 

We have shown above how the buckling of a compress¬ 
ible filament can be studied as the oscillation of a rela¬ 
tivistic nonlinear pendulum. Here are a few suggestions 
for future extensions. The instability of a flat compress¬ 
ible rod or sheet toward curved configurations is con¬ 
nected to the emergence of oscillations from rest in the 
parametric resonance of a relativistic oscillator. A change 
in the compression modulus of an elastic system (e.g., as 
a result of a temperature change) is likened to a change 
in a particle’s rest mass. Lastly, the analogy can be ex¬ 
tended to three dimensions, where the deformation of a 
compressible filament in 3D becomes analogous to a rel¬ 
ativistic rotor, as outlined in Appendix A. 


FIG. 3. Pendulum angle as a function of proper time 
(eqn (16), dashed line) and laboratory time (equations (16) 
and (19), solid line). The chosen parameters, cj )0 — 27t/3, 
co = 1, (ool/c) 2 = 0.5, correspond to a highly relativistic pen¬ 
dulum resulting in a significant time dilation. Numerical so¬ 
lution of the equation of motion (eqn (15), circles) and the 
case of a nonrelativistic pendulum (c —> 00 , dash-dotted) are 
plotted for comparison. 
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Appendix A: Analogy between compressible 3D 
Kirchhoff’s rod and relativistic rotor 


In this appendix we derive the equations of equilib¬ 
rium for a 3D compressible elastic rod under the as¬ 
sumptions of Kirchhoff’s model [18, 19]. It is shown that 
these equations are mathematically analogous to the 3D 
Euler’s equations of rigid-body rotation, where the non- 
relativistic angular momentum, L, is replaced by, yL, 7 
being the Lorentz factor. This analogy is derived while 
keeping in mind the known difficulties in the latter, rel¬ 
ativistic problem [20, 21]. 

Following the formulation in ref. [22], we consider a 
space curve described by the position vector, R(s), where 
s is the arclength parameter of the relaxed configuration. 
Defining s as the arclength parameter of the deformed 
configuration, the unit tangent vector to R(s) is given 
by, d 3 = dJA/ds. In addition, let {d^cE} be a pair 
of unit vectors perpendicular to d 3 and parallel to the 
principal axes of the filament’s cross-section. The triad 
{di, d 2 , d 3 } form a co-moving coordinate frame attached 
to the rod’s mid-axis. These vectors satisfy the relations, 


ddj 

ds 


7K X dj, 


(Al) 


where i = 1,2,3, 7 = ds/ds is the strain field, and k = 
K \ d-| + « 2 d 2 + /c 3 d 3 is the curvature vector. The local 
bending moment is given by, 


M(s) = + B 2 { 7« 2 )d 2 + B 3 ( 7K 3 )d 3 , (A2) 


force P, is given by E = / Q L e[/e(s),7(s)]ds, where 

e[/c(s),7(s)] = • (7/c) + ^(7- l) 1 2 3 4 5 6 7 8 9 +yd 3 P. (A3) 

The appearance of the 7/c term in the bending energy 
is a consequence of the requirement to keep bending 
and compression contributions independent. Minimizing 
eqn (A3) with respect to 7 and Ki (keeping in mind that 
Ki are not independent) gives [22], 

dM , . „ , 

ir-d»xp = «. m 

y( 7 -l)+d 3 -P = 0. (A5) 

In the incompressible limit eqn (A5) is redundant and 
equations (A4) become analogous to the non-relativistic 
Euler equations of a 3D rigid body, fixed at one point and 
rotating under the influence of an external force (such as 
gravity) [23, p. 200]. In this analogy, the bending moment 
takes the role of angular momentum, M -f-s> L, and the 
boundary force is analogous to an external torque, d 3 x 
P O N. Turning on compressibility effects, we have by 
eqn (A2) that M —$■ 7M. Thus, it is left to show that 
7 coincides with the Lorentz factor. First integration of 
equations (A4) and (A5) gives, 

U = ■ (7/c) - y (7 - l) 2 - yd 3 • P = const. (A6) 

Eliminating d 3 • P from eqn (A6) and substituting in 
eqn (A5) gives, 


7 = 


y/T+ ~mjy 


1 + ^/c? + ^ k \ 


B 3 „2 
Y ^3 


(A7) 


where Bi are bending rigidities. 

In correspondence with eqn (1), the energy functional 
of a 3D elastic filament, confined by a boundary constant 


This expression indeed resembles the Lorentz factor up to 
the constant prefactor, 70 = \/l + 2TL/Y , which appears 
also in the 2D problem, and which can be absorbed in 
the force P, (see discussion in the main text). 
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